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ABSTRACT 

We study models of the y-ray emission of Cyg X-3 observed by Fermi. We calculate the aver- 
age X-ray spectrum during the y-ray active periods. Then, we calculate spectra from Comp- 
ton scattering of a photon beam into a given direction by isotropic relativistic electrons with 
a power-law distribution, both based on the Klein-Nishina cross section and in the Thom- 
son limit. Applying the results to scattering of stellar blackbody radiation in the inner jet of 
Cyg X-3, we find that a low-energy break in the electron distribution at a Lorentz factor of 
~ 300-10^ is required by the shape of the observed X-ray/y-ray spectrum in order to avoid 
overproducing the observed X-ray flux. The electrons giving rise to the observed y-rays are 
efficiently cooled by Compton scattering, and the power-law index of the acceleration process 
is ^ 2.5-3. The bulk Lorentz factor of the jet and the kinetic power before the dissipation 
region depend on the fraction of the dissipation power supplied to the electrons; if it is ^ 1 /2, 
the Lorentz factor is ~ 2.5, and the kinetic power is ~ 10^^ erg s"', which represents a firm 
lower limit on the jet power, and is comparable to the bolometric luminosity of Cyg X-3. 
Most of the power supplied to the electrons is radiated. The broad band spectrum constrains 
the synchrotron and self-Compton emission from the y-ray emitting electrons, which requires 
the magnetic field to be relatively weak, with the magnetic energy density < a few times 10"^ 
of that in the electrons. The actual value of the magnetic field strength can be inferred from a 
future simultaneous measurement of the IR and y-ray fluxes. 

Key words: acceleration of particles - accretion, accretion discs - radiation mechanisms: 
non-thermal - gamma rays; theory - stars: individual: Cyg X-3 - X-rays: binaries. 



X-3 from those of confirmed black-hole binaries can be accounted 
for by Compton scattering in the very strong stellar wind from 
the compani on. That model also accounts for the lack of high 
frequencies ( lAxelsson. Larsson & Hialmarsdotteil l2009h in the 
power spectra of Cyg X-3. 

Cyg X-3 is a persistent X-ray sour ce. Its X-ray spectr a 
have been classified into five states by ISzostek et alj ( l2008h . 
who have also quantified their correlations with the radio states. 
Its y-ray emission has been discovered by the Fermi Large 
Area Telescope (LAT) and AG ILE in the soft sp ectral states 
jFenni LAT CoUaborationl l2009l hereafter FLC09; iTavani et al.l 
|2009|) . Fig. [Tl shows the average power-law y-ray emission mea- 
sured by Fermi during the active phases. The power law is rela- 
tively steep, with the photon index of F ^ 2.70 ± 0.25. Fig.[T]also 
shows the average X-ray spectrum and the 15 GHz flux during the 
periods of y-ray emission. It also shows two X-ray spectra and the 
IR fluxes in the soft state from earlier observations (see Section|2]l. 



■ 1 INTRODUCTION 

Cyg X-3 is a h igh-mass X-ray binary with a Wolf-Rayet 
(WR) companion jvan Kerkwiik et with an unusually 

short orbital period of 4.8 h, located at a distance D — 7-9 
kpc i n the Galactic plan e ( iLing. Zhang & Tan3 l2009l : iDickevI 
1 19831: [Predehl et al. liood) . Due to the lack of reliable mass 
functions and determination of the inclinati on, the nature of 
its compact object remains uncertain (see IVilhu et alj |2009| 
for a recent discussion). However, the presence of a black 
hole is favoured by considering th e X-ray and radio emission 
and the bolometric luminosity jSzostek & Zdziarskil |2003: 



and tne bolometr ic luminosity aazosteK & Z.aziarsKil l/UUa: 
Szostek. Zdziarski & McCollo ughl l2008l: iHialrnarsdotter et al.l 



20081 I2OO9I) . Also. IZdziarski. Misra & Gierlinskil ( l2010ah have 



shown that the difl^erences in the form of the X-ray spectra of Cyg 
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Figure 1. The average radio to y-ray spectrum of Cyg X-3 during the 2008 and 2009 y-ray active periods is shown in the black symbols. The error contour and 
the black dashed line show the average power-law fit with T = 2.70 ± 0.25 to the FennifLP^ data (FLC09). The black error bars show the simultaneous X-ray 
spectrum from the RXTE/ASM and 5w!/?/BAT. The black cross shows the average 15 GHz radio flux measured by the AMI and OVRO telescopes during the 
same periods. The magenta and green error bars show two other soft-state X-ray spectra, from BeppoSAX and RXTE, respectively. The magenta dashed line 
shows a 1.5-100 GHz, Fy oc v""-', spectrum, similar to the spectra measured in that range during a 2001 radio outburst of Cyg X-3, and the two magenta 
crosses show IR measurements in a radio-flaring state. See Section|2]for details. 



We see that the y-ray power law spectrum has to have a low- 
energy break somewhere between ~1 and ~100 MeV in order 
not to produce more X -ray emission than observed. Furthermore, 
IZdziarski et alj ^201 j) show that the hard X-rays up to at least 
100 keV during the intervals with y-ray emission have the orbital 
modulation pattern characteristic to wind absorption and scatter- 
ing. On the other hand, the GeV orbital modulation is shifted in 
phase with respect to the X-rays by ~ n/2 (FLC09). This appears 
to imply that the contribution of the spectral component observed 
in the GeV range to the 100 keV flux is at most weak. The GeV 
power law emission and its orbital modulation appear to be due to 
Compton up-scattering of the stellar emission from the companion 
WR star by relativistic electro ns in the jet of this source (FLC09; 
iDubus. Cerutti & Henrill2010bl . hereafter DCHIO). 

In this work, we first calculate the average X-ray and radio 
emission during the y-ray active periods. Then, we study emission 
due to Compton up-scattering of blackbody photons by power-law 
electrons with a low-energy break. Finally, we apply our theoretical 
results to the broad-band X-ray/y-ray (Xy) spectra of Cyg X-3, and 
obtain strong constraints on the electron distribution in the y-ray 
emitting region and on the parameters of the jet. 



2 THE RADIO-X-RAY SPECTRA 

We use X-r ay monitoring data from the Swift Burst Aler t Tele- 
scope (BAT: iBarthelmv et alj[2003 : iMarkwardt et al.ll2005h in th e 
form of a 14—195 keV 8-channel light curve ( Zdziarski et alj2012h . 
The channels are between energies of 14, 20, 24, 35, 50, 75, 
100, 150 and 195 keV. Also, we use X-ray monitoring data from 

i3. 



the the All-Sky M onitor (ASM;,Bradt, Rothschild & Swank 1993' 
lLevineetalJl996h on board Rossi X-ray Timing Explorer (RXTE). 
The ASM has three channels at energies of 1.5-3 keV, 3-5 keV 
and 5-12 keV. In the radio range, we use the 15 GHz data from 



the AMI Large Array and the Owens Valley Radio Observatory 
(OVRO), which data were used in FLC09. 

We then calculate the average fluxes in each of the X-ray 
ch aimels (converting the count rates to fluxes using the method 
of Izdziarski. Poolev & Skinneill201ll) and at 15 GHz during the 
two y-ray active periods for which the average y-ray spectrum 
of FLC09 was obtained, which are MJD 54750-54820 (in 2008) 
and MJD 54990-55045 (in 2009). The resulting average X-ray 
spectrum and the average radio flux simultaneous with the y-ray 
spectrum are shown in Fig. [T] The average 15 GHz flux equals 
0.38 ± 0.04 Jy, with a large rms of 0.43 Jy, reflecting a strong flux 
variability from a few mjy to ~2 Jy. 

To illustrate the likely form of the radio spectrum during the 
y-ray active periods, we plot a radio power law between 1.5 GHz 
and 100 GHz with an energy index of 0.5 (characteristic to un- 
cooled optically-thin synchrotron emission), which is in the middle 
of the ~0.4-0.6 index range measured i n that frequency range dur - 
ing a 2001 radio outburst of Cyg X-3 bv lMiller- Jones et al.l |2004|) . 
To show the likely form of the broad-band soft-state radio-flaring 
spectrum of Cyg X-3, we also show two IR measurements, at 4.5 
yum and 11.5 yum, taken in the flaring radio state on 1997 June 18, 
with the quiescent-state fluxes (w hich are, most lik ely, due to the 
stellar wind emission) subtracted dOglev et alj|200lh . The 15 GHz 
flux measured by the Ryle telescope at the time of the IR measure- 
ment, MJD 506217.(39-42), was 0.64 Jy, i.e., about 1.5 the average 
value during the y-ray active periods, and the 8.3 GHz flux from the 
Green Bank Interferometer was 0.59-0.72 Jy. ( Note that the MJD 
of the measurements given in lOglev et al.ll200ll are in error.) 

In Fig. [T] we also show one of the average soft-state 
X-ray spectra from RXTE Proportional Counter Array (PCA) 
and the High Energ y X-Ray Transient Experiment (HEXTE) of 
ISzostek et alj ( 1200 8l) . nam ely their spectrum #4, and th e soft spec- 
trum from BeppoSAX of ISzostek & Zdziarskil ( |20C)8|), which is 
similar, but of better quality, to the spectrum #5 of Szostek et al.l 
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( l2008h . We see that our average ASM/BAT spectrum lies between 
these two spect ra, thus it is in termediate between the soft and ul- 
trasoft states of lSzostek et alj ([2OO8). As shown bv lZdziarski et alj 
there is a close correspondence between the canonical X- 
ray states of black-hole binaries and the X-ray states of Cyg X-3. 
The soft and ultrasoft states in black-hole binaries are dominated by 
emission of an optically-thick accretion disc up to an energy of sev- 
eral keV, and by a power-law like tail, probably of coronal origin, at 
higher energies. The same situation is then, most likely, present in 
those states of Cyg X-3, with a modification due to the passing of 
that emission through the very strong stellar wind of the donor WR 
star. Some jet contribution is also possible at hard X-rays, which 
issue is studied in this work. 



3 ANISOTROPIC COMPTON SCATTERING 

In this Section, we consider Compton scattering of soft photons by 
a cloud of relativistic electrons isotropic in the jet comoving frame. 
Thus, the photon energy, electron Lorentz factor and the scattering 
angle below are given in this frame. 

The problem of Compton scattering of a mono-directional 
photon beam by a cloud of relativistic electrons with a Lorentz 
factor 7 » 1 and an i sotropic angular distribution into a given an- 
gle has been solved bv lAharonian & A tovan ( 1981). Their equation 
(20), valid from the Thomson to the Klein-Nishina regimes, gives 
the flux per electron and per solid angle, which can be written as. 
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where 6o and 6 are the energy of the incoming and scattered photon, 
respectively, in units of , i? is the scattering angle, ho [cnT^ s" ' ] 
is the number flux of incoming photons, is the electron mass, 
and CTj is the Thomson cross section. Also, 



60 « 6 < 



(3) 



1 -I- 2x£oJ 

where the former constraint expresses the applicability of the » 
1 condition, and the latter, equivalent to r < 1, is kinematic. 

The above rate can be then integrated over an electron dis- 
tribution. We consider the case of a power-law distribution with 
cut-offs. 
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(4) 



where the constant K specifies either the electron density or their 
total number. Accounting for r < 1, we have. 
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Figure 2. An example spectrum from a cloud of electrons irradiated by a 
beam of monoenergetic photons with eo = 5 X 10"^, no = 1, emitted at 
?? = 90° with respect to the beam for p = 4, yi = lO', 72 — > °°, and 
the normalization corresponding to K = 1 . The red solid and green dashed 
curves correspond to the Klein-Nishina formula (6) and the Thomson limit 
with a sharp Klein-Nishina cut-ofF, equations )A1HA4K respectively. 



where B,, is the incomplete beta function. In terms of e, the first con- 
dition in equation Q reads 6 < 2xeo7[/(l + 2x6o7i). Equation ^ 
for the second case in equati on Q is equivalent to equation (33) of 
lAh aronian & Atovad i fl98ll) . For integer or half-integer j, By( j, 0) 
can be expressed relatively simply by elementary functions, and, in 
general. 



y<l. 



(8) 



For a finite 72, the flux can be obtained by subtracting the rate with 
72 substituted for 71 in equation Q from the rate of equation 
Then, the spectrum will be null for e > 2xeoy\l{\ + 2xeoy2). Fig. 
|2] shows an example of the spectrum for parameters roughly appli- 
cable to Cyg X-3, and compares it to the Thomson-limit spectrum 
of Appendix|A] We see that for the chosen parameters, the average 
slope of the actual spectrum for a decade above the break (corre- 
sponding to 7i) is substantially steeper than of that in the Thomson 
limit. 

Equation l|6j can be then integrated piece-wise over a distribu- 
tion of irradiating photons, ho(eo). 
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where q <g: 1 is a constant assuring the eo ^ f condition, at 
which equation (TJ is valid. When K corresponds to electron den- 
sity, ed;j/dedn is the emissivity. When K corresponds to the total 
number of electrons, it gives the total photon production rate, and 
the dimension of equation l|9) is s"' . In the latter case, the differen- 
tial luminosity and the observed flux (neglecting relativistic correc- 
tions) are given by. 



dL 



edh 
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1 dL 
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(10) 



dedn dedn 
respectively. 

The seed photons can be, in particular, from a blackbody emit- 
ter (e.g., a star) with a radius, R,, and temperature, T. When the 
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Figure 3. The thick sohd bkie, green and red curves show example Klein- 
Nishina spectra in the cloud frame from a stationary (jSj = 0) cloud of elec- 
trons irradiated by a beam of blackbody photons with the emission at an an- 
gle with respect to the beam given by cos & = -0.8, 0, and 0.8, respectively, 
for /) = 4, yi = 10^ ji ^ oa,T = YO^ K = i,R = R.. The correspond- 
ing dashed curves show the Thomson-limit spectra, see Appendix |A] The 
black dotted line shows the spectrum of equation jA9) for cos i? = 0.8. 



distance of the electron cloud from the stellar centre, R, is » R,, 
the blackbody photons form an almost mono-directional beam, in- 
cident on the electrons. Here, we take into account that the electron 
cloud may be located in a jet or counterjet moving with respect to 
the star, for which the Doppler factor is. 



1 



rj(l-/3je.-ej) 



(11) 



where /Jj is the jet velocity, Fj = 1/(1 - ySj ) is the jet Lorentz 
factor, and e, and ej, are the unit vectors along the direction from 
the star towards the electron cloud, and along the jet, respec- 
tively, see fig. 1 in DCHIO. Then, the soft photon energy in the 
stellar {= observer) frame is e^I),, and this energy has a black- 
body distribution at kT/m^c ^. Since ^!o(eo) is a relativistic invariant 
( iBlumenthal & Gould|[T97ol) . 



2n 
7^ 



(I) 



(12) 



where k and h are the Boltzmann and Planck constants, respec- 
tively. 

Fig. [3] shows example spectra obtained using equation l[9} 
for blackbody irradiation, also comparing it to the corresponding 
Thomson-limit spectra of Appendix IaI Such a model can give the 
low-energy break required by the broad-band Xy spectrum ob- 
served from Cyg X-3, see Fig. [T] We see that the photon break 
energy moves to lower energies as i? decreases, as implied by equa- 
tion 0. The low-energy parts of the spectra are almost independent 
of the angle. Similarly to the case of mono-energetic incident pho- 
tons, Fig. [21 the Thomson-limit spectra above the break have the 
slope significantly harder than the actual spectra, but the difference 
between the two decreases with the decreasing scattering angle, &. 



4 OBSERVED SPECTRA 

Here, we calculate the spectra in the observer's frame. We consider 
a steady-state jet, in which the observed emission comes from a 
given spatial range in the observer's frame. This is compatible with 
the dynamical time scale of the jet of the order of tens of seconds 
(see Section 15. 2t whereas the y-ray emission was detected over 
time scales of days/weeks. Furthermore, the strong orbital modu- 
lation in y-rays (FLC09) requires the y-ray emitting region to be 
approximately stationary. 

The jet and counterjet Doppler factors with respect to the ob- 
server are 



1 



Fj(l -/3jeob.s-ej) 



1 



rj(l H-^jeobs-Cj) 



(13) 



respectively, where Bo bs is the unit vec tor towards the observer We 
use the result of Dubus. Cerutti & Henri t2010a,). 



1 - cos J? = DjD. (1 - Cobs-e.) , 



(14) 



where tabs • e, is the cosine of the orbital-phase dependent angle 
between the direction from the star to the jet and from the jet to the 
observer in the observer's frame (which angle in the jet frame is the 
scattering angle, ff). For the jet emission. 



dF 
dE 



edn 



DTj dedfi 



(15) 



where E is the observed dimensional photon energy, and the 
flux transformation to the observed frame is for a steady-state jet 
dSikora etal.ll997h . We note that the form of dF/dE above assumes 
the the energy units in F and E are the same, whereas they are of- 
ten assumed to be different (e.g, erg and eV), which, however, can 
be easily accounted for For the counterjet, Dq and the correspond- 
ing D,, R should be used, and the two observed fluxes should be 
added. Given the observed spectrum, this transformation yields the 
normalization of the emitting electron distribution corresponding to 
the actual number of electrons in the considered jet region. On the 
other hand, the transformation used in DCHIO, with instead of 
Dj/Tj, corresponds either to emission of a single moving blob or 
to the observed number of electrons in the emitting part of the jet. 
Note that in the latter case the observed number is different in the 
jet and the counterjet, which was not accounted for in the treatment 
used by DCHIO. 

For a pure power-law emission from Thomson-limit scatter- 
ing, as in equation l |A9t . applying the above relations results in the 
energy flux as given by equations (1) and (3) of DCHIO (except for 
the form of the dependence on Dj), with no dependence on O,. In 
a general case, we need to apply formulae ( ll4H15t to equation ^ 
(i.e., substitute e and x as above and multiply the rate by D^/Tj), 
integrate it numerically, and then repeat the procedure for the coun- 
terjet. 

The distance between the electron cloud and the stellar cen- 
tre, R, and the components of the vectors e,, eob.s and Cj are given 
by equations (4) and (5) of DCHIO, respectively. They depend on 
several parameters of the system, namely, the orbital separation, d, 
the height of the electron cloud along the jet, H, the binary incli- 
nation, (, the inclination of the jet with respect to the normal to the 
binary plane, given by the azimuth, dj and the polar angle, tpj, and 
the orbital phase, 9 (note that DCHIO use a non-standard definition 
of 0, with 6 = ±n/2 rather than the usual 0, n, at the conjunctions). 
We include here both the jet and the counterjet, for which the polar 
angle is (f>j+n (and its unit vector is -ej). We note that the model of 
DCHIO neglects eclipses of the counterjet by the star. We neglect 
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Figure 4. Example spectra for the best-fit black-hole model of DCHIO in 
the observer's frame, but using the full Klein-Nishina spectra, p = 4, and 
modified to account for the flux transformation of equation (15), see Section 
|4] The dashed red and blue curves give the spectra from the jet at the orbital 
phases (defined as in DCHIO) of = Q.2n and 0.7;r, respectively, and the 
corresponding dotted curves give the spectra from the counterjet. The solid 
black curve gives the spectrum from both the jet and counterjet averaged 
over the orbit. The parameters ai'e y\ = 10-', y2 ^ oo^ T = 10^ K. The 
normalization corresponds lo K = 1; then the unit of vertical axis is keV 



them here as well for consistency with the adopted assumption that 
both the blackbody emission and the scattered emission are point- 
like. 

We use here the model of DCHIO with their assumed black- 
hole binary parameters, with = 4.1 x 10" cm, / = 30°, T = 10' 
K and R, = 2.3Rq. However, since we use a different relativistic 
transformation between the jet and observer's frames, we have re- 
fitted the model to the observed orbital modulation of y-rays using 
equation (Tsj. We have obtained // = 7.6 x 10" cm, dj = 319°, 
0j = 39°, which are the same as in DCHIO, but /3j = 0.47 
(Fj ^ 1.13), which is somewhat higher than their value of 0.41. 

At these parameters, relativistic beaming is moderate, Dj ^ 
1.50, ^ 0.62,0. ^ 1.45-1.60, l-eobs*. ^ 0-0.4 for the jet and 
1.6-2 for the counterjet, and R/R, ^ 3.6-6.7. Since we use here the 
Klein-Nishina cross section instead of the Thomson approximation, 
the actual spectrum is softer than in the Thomson limit, where p = 
2r - 1 = 4.4. Here, we use p = 4, which approximately reproduces 
the observational best-fit spectrum. 

Fig.|4]shows two example spectra for the electron distribution 
with yi = 10^^ and d = 0.2/r and 0.7;r. We see they look similar 
to those in Fig. [3] except for an additional shift of their relative 
normalization, introduced by the phase-dependent D, and R/R,. 
We then average the observed spectrum over the orbital phase, 9, 
which we plot in the solid curve. 

We note, however, that we have to take into account the elec- 
trons below yi. Even if the electrons are accelerated only above yi, 
they lose energy via Compton, synchrotron, and adiabatic losses 
and form a distribution below yi. Hereafter a dot will denote a time 
derivative in the jet frame (as in Section[3]l, and d.../dt will denote 
a time derivative in the observer's frame. The loss rates in the jet 



frame are given by 

4,/kNO"Tf^radr" 
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(17) 



respectively, and is the blackbody energy density within the 
electron cloud, /kn < 1 gives the Klein-Nishina reduction with re- 
spect to the Thomson limit, B is the magnetic field strength, and the 
factor of 2/3 in y^d accounts for the expansion being in two dimen- 
sions only. The cooling rate for a single electron for monoenergetic 
se ed photons us ing the Klein-Nishina cross section was calculated 
bv lJonesl ( [T96^,ri968il . The Lorentz factor at which the Compton 
and adiabatic rates equal each other, and the electron distribution 
has a break, is 
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o-jft^ikTY H 
which for the assumed parameters equals to, 
T 0.018 
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(io'k) 
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where /kn = 1 was assumed and the factor of 0.018 is the value of 
the orbital overage of (R,/RD,)^ in f/pad- Thus, the break Lorentz 
factor at our parameters is well below the minimum Lorentz factor 
of ~ 10^ required to explain the observed y-ray spectrum. We note 
that the dependence of yt, on H and y6j is rather complex; /f"' ap- 
pears in equation J18t . but R also depends on H and D, depends on 
H and /Jj. Furthermore, the parameters are mutually connected via 
the requirement of fitting the observed orbital modulation. 

We assume that the electrons are accelerated at a power-law 
rate, which in either the jet or counterjet frame is given by. 



yi < y < j2- 



(20) 



where A^i„j is the normalization factor. Hereafter, K and A"i„j corre- 
spond to the total number of electrons in the jet (and not to their 
density). Then, assuming Compton losses in the Thomson limit 
(/kn = 1) and synchrotron losses, and for yi > yj, (fast cooling), 
the steady-state distribution (in either jet or counterjet) will approx- 
imately be. 



Ky-i\ 
0, 



(21) 



yo < y < jh, 
yb < y < yi; 
yi < y < y2; 
otherwise, 

where yo ~ 1 is a minimum overall Lorentz factor. The steady- 
state electron kinetic equation in the comoving frame is N{y) = 
y"' 2(y)dy, where y is total loss rate. Assuming the dominant 
losses above yi are Compton in the Thomson regime, A^i„j is then 
related to K by, 

4o-jU,,i{p-2)K 



(22) 



Fig. |5ja) compares the models with p = A with the obser- 
vations. Here, we impose K to match the best-fit Fermi spectrum 
and require that both the low-energy break is at ^0.1 GeV, and the 
X-rays at ~100 keV are not overproduced. We show the average 
spectra for the electron distribution of equation M\\ with yb = 130 
and yi = 700 and 1500. As expected, the low-energy break is at an 
energy oc yj. Fig. |5jb) compares the models with p = 3.5 and the 
normalization somewhat below the best fit, which both are within 
the observational uncertainties. We see that these results imply a 
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Figure 5. The radio to y-ray spectra spectra (the same as in Fig.[T) compared to the Compton and synchrotron models for the electron distribution of equation 
(21) with 7b = 130, 72 — > oo, r = 10^ K. (a) The red dashed and blue solid curves show the Compton spectra for yi = 700 and 1500, respectively, and p = A. 
The model spectra are normalized to the best-fit average Fermi spectrum. The dashed blue curve shows the synchrotron spectrum from the y-ray producing 
electrons for yi = 1500 and the magnetic field given hyi]B = 10"^ The blue dotted curve shows the synchrotron self-Compton spectrum, (b) The same for 
p = 3.5 and yi = 300 and 1300, for the red dashed and blue solid, dashed, and dotted curves, respectively. The normalization of the solid curves is below the 
best fit but within the error contour. We see that the value of the electron low-energy break is constrained to 300 <yi < 10^, but models with yi < 10^ give 
strong contribution to the hard X-rays. The cyan dotted curve shows the Compton spectrum corresponding to the electron distribution of equation (23) with 
71, = 10^ and yi = 10. This model is in principle possible, but the value of yb used for this model is several times our estimate, equation )19t . 



constraint of 300 ^ 7i ^ 1500. However, models with 71 < 10^^ give 
strong contribution to the hard X-rays, which appears in conflict 
with the related result that the orbital modulation at ~ 100 keV 
during the y-ray emitting intervals is characteristic to bound-free 
absorption and Compton scattering by the stellar wind and o ut of 
phase with the > 0. 1 GeV modulation jZdziarski et alj|2012h . We 
note that the index of the accelerated electrons in our models is 
p-l^ 2.5-3. 

The red and blue curves in Fig. [6] show the orbital modulation 
pattern at 100 keV and 0.2 GeV of the model shown Fig.[5ja) with 
p = 4 and 71 = 700 (at 0.2 GeV, the modulation pattern of the mod- 
els with 71 = 700 and 1500 are identical). We see that the patterns 



at the two energies are very similar, and have the maxima around 
the superior conjunction, whereas the X-rays have the minimum at 
this phase. Given that the optical depth through the wind from the 
electron cloud is much lower than that from around the compact 
object, we can use the orbital modulation of X-rays to distinguish 
the X-ray source location close to the compact object from that in 
the scattering cloud in the jet. 

On the other hand, the parameters used in equation ([19), e.g., 
T, R,, bear large uncertainties, and we cannot exclude in principle 
a much larger value of 7b. For 71 < 7b (slow cooling), we would 
then have approximately, 
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Figure 6. The orbital modulation of the (fast-cooling) model with yi = 700, 
yb = 130, and p = 4 (which spectrum is shown by the red dashed curve in 
Fig. [5^) at 10" keV (red solid curve) and 0.2 GeV (blue solid curve). The 
blue dashed and dotted curves show separately the contributions of the jet 
and counterjet at 0.2 GeV. The modulation pattern at 10- keV of the (slow- 
cooling) model shown by the cyan dotted curve in Fig. |3Ib) is relatively 
similar to that shown by the blue solid curve. The superior conjunction 
(compact object behind the WR star) is at B/ln = 0.25, around which the 
X-rays show the minimum due the maximum absorption/scattering. On the 
other hand, our models have the maxima around this phase, which appears 
to rule out any substantial contribution of the jet to the X-rays. 



Niy): 



Kj-i', 
0, 



i-p 



Jo <j <yi 
7i < 7 < 7b 

7b < 7 ^ 72 
otherwise, 



(23) 



The cyan curve in Fig.[5jb) shows the case for p = 3.5, 7b = 10', 
7i = 10. We consider this model relatively unlikely given the es- 
timate of 7b of equation l |19l l, and consider below models with 
7b < 7i. Although it reproduces well the high-energy tail of the 
softest X-ray spectrum, its orbital modulation is very similar to 
those shown in Fi g. [6] which is in confl ict with the observed hard 
X-ray modulation dZdziarski et al.ll2oT^ . 



5 THE JET STRUCTURE 
5.1 The jet power 

The amplitude of the 7-ray orbital modulation close to unity indi- 
cates the 7-ray emission region is rather compact, and the jet does 
not emit along a range of heights large compared to the average 
distance of the source from the compact object, H. Otherwise re- 
gions at different heights would have different modulation patterns, 
strongly reducing the net modulation amplitude. Thus, prior to the 
7-ray emission, we assume non-radiating electron-ion jet and coun- 
terjet moving with the Lorentz factor of Fj q. The sum power of the 
jet and counterjet is then dominated by the bulk motion of cold 
ions. 



Pj.o = 2m,c2(Fj.o - 1) 



dM 
"dT' 



(24) 



where dNJdt is the ion number flux in either the jet or coun- 
terjet in the observer's frame and nij is the ion mass. Given 
that Cyg X-3 contains an He donor jvan Kerkwiik et all 1 19961 : 
iFender. Hanson & PoolevI 19991) . nii ^ 4mp, where is the proton 
mass. The jet then enters a shock region, reducing its bulk Lorentz 



factor to Fj, which for the best-fit parameters of the orbital mod- 
ulation model is Fj ^ 1.13 (Section|4ll. The fraction of the initial 
energy dissipated is, 

rj,o-Fj 



1 j.o - 1 

The power supplied to the electrons is given by. 



(25) 



(26) 



where r]^ is the fraction of the dissipated energy supplied to 
electrons in the shock acceleration region. We assume (e.g., 
ISpitkovskvll2008l) that all electrons in the dissipative zone are ac- 
celerated/heated to relativistic energies. Thus, 

dM _ w, diVe dTVe 
dt n, dt ' df 



r.l 



(27) 



where dN^/dt is the total electron number flux in the observer's 
frame, Q{~y) is the electron acceleration/heating rate in the jet frame, 
and the ion/electron density ratio is iijn^ ^1/2 for He in the ab- 
sence of pair production. 

Combining equations | |24H27| |. we find, 



rj,o = Fj 1 + 



m^n^ 7in 



(28) 



where 7i„j is the average Lorentz factor of the accelerated/heated 
electrons in the jet frame. Here, we approximate Q(y) , consisting 
of a r elativistic Maxwellian and a power-law tail (e.g.. lSpitkovskvl 
l2008h . as a power law with the index of p - 1 and a low energy 
cut-ofi', see equation l l20t . Then, for 72 — » 00, 



/ e(7)ydy ^ 71 (p- 2) 
/ e(r)d7 " P - 3 ' 



dWe 

"dT 



Fj(p-2)' 



(29) 



for p > 3 and p > 2, respectively, and where A^i„j is given by 
equation l l22b . Equation l l28t then yields, for p = 4, Fj,o Fj -1- 
(0.62/j/e)(7i/10'). The value of t]^ remains unknown; if the dissi- 
pated power is divided equally among electrons and ions, r\^~ \ 12. 
For this value and 71 = 10', Fj o - 2.38 and rji^si - 0.90. 

The radiated power, Piad, is related to the power in the elec- 
trons by, 



^niP c. 



(30) 



where zymj is the radiation efficiency. We then calculate the initial 
jet power and the radiative power as 



Fj - 1 m,n, 71 (p - 2) 



Fj m^n, 



Sj]-"Ka-jcU,,i 



Pnid - '/rad- 



(P - 3)77, 

(p - 2) SK<ncUr,6 



P- 



(31) 



(32) 



Note that in the above derivation we did not need to specify the 
extend of the dissipation zone. Since 7b <k 71 , see equation l |19l l, 
the radiative efficiency in the dissipation zone is 7/,ad - 1. 

The normalization constant in Fig. \5la), for p = 4, is 
K/D- = 2.1 X 10\ which implies K ^ 1.0 x 10^°(Z)/7 kpc)^ 
(For the models in Fig. 0) with p = 3.5, K/D^ = 2.2 x 
10'.) This yields the powers averaged over the orbital phase of 
P,o - 7.2[0.11 -K0.5/;7e)(7i/10')](7i/10')--''(£»/7kpc)^10" erg 
s- ' , and Pi-ad 3 .6(71 / 1 0' )'-''(D/7 kpc)- 1 0" erg s- ' . (We note that 
the jet radiative output of electrons with 7 > 10' given in DCHIO 
is mistakenly too large by a factor of 4n.) The Klein-Nishina cor- 
rected radiative power for the steady-state electron distribution of 
equation l l21t is slightly lower. 



© 2012 RAS, MNRAS OOO.fTirB] 



8 A. A. Zdziarski et al. 



We can compare the total mass flow rate in the pre-shock 
jet, i.e., including the rest mass and the associated kinetic en- 
ergy, Mj = 2mjrjodA'j/df, to the mass accretion rate estimated 
from the bolometric luminosity of the source. At rj^ = 0.5, 
Mj 1.5(O/7kpc)H0" g s"'. On the other hand, the bolo- 
metric luminosity of Cyg X-3 in the soft state calculated by 
ISzostek & Zdziarskil Jiooj) is ^ 2(D/7 kpc)^10^'* erg s"', which, 
at an accretion efiiciency of eaccr = 0.1, requires Maccr - 
2(£)/7kpc)^10'* g s"'. Thus, a relatively small fraction of the 
mass flow rate at the outer boundary of the accretion source in 
its soft state is sufficient to power the jet, unless rj^ ^ 1. We 
note that this fraction is still much higher than that estimated 
for the hard- state jet in the black- hole binary XTE Jll 18-480 
of ~ 0.01 bv lYuan. Cui & NaravanI l2005h . The Eddington limit 
on the mass accretion rate for He corresponds to Me - 3 x 
10''(M/10Mo)(6,,,,/0.1)-' erg s"', where M is the black-hole 
mass. Adopting this limit imposes a constraint on r;^. 



5.2 Electrons, cooling and pair production 

Based on the steady-state electron distribution of equation l l21t . we 
can calculate the total number of relativistic electrons and their total 
energy in the comoving frame for the distribution of equation i2l\ 
(for 72 — * °°) for either the jet and counterjet. 



1 + In — 

70 



P- 



P- 



P>U 



P- 



1 , yi 

— +ln — 

2 7b 



p>2. 



(33) 



(34) 



respectively. This, for p - 4, 71, = 130, 71 = 10', yields ^ 
4.4(£)/7kpc)M0^^ and ^ 2.1(O/7kpc)-10'* erg. Then, the ra- 
diative time scale averaged over the electron distribution in the jet 
frame is 



1 103/ llO'K/ 



(35) 



The cooling time of an individual electron, ~ (lO^'/y) s [neglecting 
/kn. see equation fT6l>1. is approximately equal to (ty^i) at 7 = 71, 
which is close to the average electron energy, see equation J29l >. 
Also, /3jrjc(fiad) also gives the minimum size of the emission re- 
gion. It cannot be more compact because the electrons would not 
then have time to cool. 

The dynamical time scale, in the jet frame, is longer (which 
also follows from 7b < 71 for our parameters), 

H I H 



fdyi 



50 



( ] 

\8x 10" cm/ 



(36) 



Since the emission at > 0.1 GeV is due to electrons with 7 > 10', 
their cooling time being < t^y^ is compatible with the orbital mod- 
ulation of photons with energies > 0. 1 GeV being close to 100 per 
cent, requiring the corresponding emitting region to be compact. 
On the other hand, low-energy electrons will radiate and lose their 
energy over longer ranges of the jet length, and thus the depth of 
the orbital modulation is expected to decrease somewhat with the 
decreasing photo n energy at £ < 0.1 GeV. 

As shown bv lCerutti et alj I2OI ih . pair production on accre- 
tion blackbody disc photons in the 7-ray emission region is negli- 
gible. This process would absorb only 7-rays emitted from a vicin- 
ity of the accretion disc, at distances < 10^ cm for the parame- 
ters adopted here, or, at < 10'" cm if the disc emission get full y 
isotropized by the stellar wind (but see lSitarek & Bednarekll2012h . 




10= 
E [GeV] 

Figure 7. The optical depth to pair production on stellar photons for y-rays 
produced in the electron cloud. The red solid and blue dashed curves are for 
6 = 0.9 and 2.0 (in radians), respectively. The optical depth is around the 
maximum for the former phase, and it goes to a minimum of Tyy « 1 at all 
energies for 9 ^ 3.2. The phase of 9 = 2.0 is intermediate. The distance from 
the stellar centre and the angle to the observer with respect to that direction 
are R = 8.5 X lO' ' cm and 6.2 X lO" cm, and arccos(e,-eobs) = 54° and 40° 
at 6 = 0.9 and 2.0, respectively. 



The present model does not predict 7-ray emission in these regions. 
At higher energies, 7-rays at > 10 GeV are above the threshold 
for pair production on stellar photons. The degree of attenuation 
strongly depends on the assumed inclination (which, in turn, is re- 
lated to the fitted value of H since the observed modulation depth 
close to unity requires that the star, the electron cloud and the ob- 
server are aligned at some phase.) For the parameters used here, the 
maximum optical depth, Tyy, to this process is moderate, as shown 
in Fig. [7] which has been calculated using the method of iDubuj 
( l2006al) (taking into account the finite size of the star). However, 
an inclination > 30° would also yield a lower value of H, with both 
changes significantly increasing Tyy. 



5.3 Magnetic field 

We then consider the magnetic field in the 7-ray emitting region. 
We assume that the magnetic energy flux in the downstream region 
is a fraction, 77^ < 1, of the dissipated power. Including both the 
jet and counterjet, the magnetic field in the jet comoving frame is 
given by, 

R- 

— - ■■■■ - (37) 



'7iad'7c 



where Oj is the jet opening angle in the dissipation region, 77„d - 1, 
and the magnetic energy flux in the second equality is expressed in 
terms of the quantity closest to the observations, i.e., Praj. 



For Prad estimated in Section O B l20{r]B/n^y'^&^' G, 
and B-/8n ^ 56O(77g/77e)0j"" erg cm"'. In comparison, the black- 
body energy density within the electron cloud is t/,ad - 3.4 x 
10^(T/I(f K)'* erg cm"\ see equation l ll6t . The magnetic field is 
constrained by the contribution of the synchrotron component to 
the broad-band spectra, which results in an upper limit on tjb, dis- 
cussed below. 

The opening angle of the jet of Cyg X-3 on a 10 mas 
scale [(ti/7kpc)10 ' ^ cm] based on 2001 radio-outburst data of 
iMiller- Jones et in l l2004h has been estimated as 0j = 5.0 ± 
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0.5° by iMiller- Jones Fender & Nakail ( |2006|) . We note that 
iMiller- Jones et alj i 2004h have fitted the jet curvature they observed 
as due to a jet precession with the precession angle of 2.4°, which 
precession may smear the observed image, with the actual opening 
angle possibly being < 5°. 

We approximate here the synchrotron spectrum using a delta- 
function approximation, in which the energy of a synchrotron pho- 
ton averaged over the pitch angle, a, is given by, 



1, B„ 



eh 



(38) 



where 5^ is the critical magnetic field. This formula follows from 
the co rrespondence between the synchrotron and Compton pro- 
cesses telumenthal & Gou ld' 197^. in which synchrotron emission 
is considered to be Compton scattering of virtual photons at the di- 
mensionless energy of 5/5cr- We then require that the power of the 
synchrotron emission in this approximation equals the actual syn- 
chrotron power, which then yields, 



6dns o-jBi(B/B,,y 



dedQ. 



48;r2;i 



-N 



B 



(39) 



Since N{y) gives the volume-integrated electron distribution, this 
formula gives the synchrotron emission from the entire source, 
analogously to the treatment in Section|3] We have found that this 
formula provides a good approximation to the actual synchrotron 
spectrum. For a power-law spectrum, N{y) = Ky^'', with a value 
of p > 1/3, the spectrum averaged over the pitch angle is (cf. 
Ijones. O'Dell & SteirJI 19741) , 



edns 



■ C 



Ci = 



r(2£-i)r(^)r(i^) 



(40) 



(41) 



where F is the gamma function. In the delta-function approxima- 
tion, equation l |39l l, C\ is set to 1. This approximation is fully accu- 
rate for p = 3, for which Ci(p) = 1. For p = 2 and 4, C[(p) = 1.14, 
1.20, respectively. For a broken electron power law, the break en- 
ergy also appears close to the actual one. We note that a 1, which 
provides a phenomenological best fit to accurate results, does not 
correspond to the average synchrotron emission of a single elec- 
tron, for which a = 2^3"''^5"' sincf. 

We need to take into account self-absorption of the syn- 
chrotron spectra. The synchrotron self-absorption coefficient in the 
jet frame for an electron power-la w distribution aver aged over the 
pitch angle can be expressed as (cf. |jones~et alJl974h , 



ttS = C; 



TTCTj K 



— -y Nir) cm , 



C2 = 



3^r(^)r(2^)F(^) 



(42) 



(43) 



where Of is the fine-structure constant, V is the source volume in 
the jet or counterjet frame, C2 - 1 for p = 3. The second formula in 
equation i42l gives ag in the monochromatic approximation, with 
y(e) given by equation OSt . This also gives as corresponding to 
emission by electrons with a given 7. The volume depends on AH, 
the length of the emission region along the jet in the observer's 
frame, and on the jet radius, @jH. The value of AH is relatively 
uncertain; it is >y6jFjc(frad) - 0.24// based on the cooling argument 
(Section [5.2t . and AH « // to account for the depth of the orbital 



modulation close to unity (DCHIO). Given these constraints, we 
adopt 



AH = ySjFjcXf^d); V ^ nQ-H-AHT,. 
For our adopted parameters, V 2.8 x 10^'' 



(44) 



cm . For B given 
by equation l l37t . we can determine the turnover energy (in the 
jet frame), at which the optical depth through the jet spine in 
the observer's direction, Ts(ei) = 2c!'s0j///2)j sin; = 1. We 
find this takes place between the radio and IR ranges, around 
~ 0.007(0.5;7B/10"^77e)''"' eV, and in the part of the spectrum emit- 
ted by electrons dominated by adiabatic losses, N{y) oc y"' of equa- 
tion i2l\ . Below El, «s oc e^'^. 

The resulting spectra, from both the jet and the counterjet, 
and taking into account the relativistic transformation of equation 
l |15l l, are shown in Figs. [Sja-b) for rjg = 10"^. The synchrotron 
spectrum has the shape similar to that of the Compton one, and, 
for the best-fit parameters, its peak, from electrons with 71, is at 
6s - 16(7;B/77e)'''(7i/10')- eV. We see that IR measurements si- 
multaneous with those in 7-rays would provide constraints and/or 
a measurement of the jet magnetization, r]B- If the shown IR mea- 
surements during radio flares are representative for 7-ray active pe- 
riods, the jet is relatively weakly magnetized, as an increase of B 
would increase the synchrotron flux at the peak oc g<''+i)/2 jjQjg 
that r]B ~ 10"' is consistent with the t heore ti cal estimates for m ag- 
netized shocks oflMedvedev & Loebl ( ll999l) . lMedyedevl ( I2OO6I) and 
ISironi&Spitkovskvl ( l201 l|). At the above parameters and 71 = 10 , 
/J ^ 60 G. As a consequence of AH derived from electron cooling, 
the ratio of the magnetic field energy density to that in the electrons 
equals to rjB/rj^, i.e., (B^/8n)/(Ec/V) = r/s/rj^. The field strength is 
thus much below equipartition. 

We also need to consider the synchrotron self-Compton pro- 
cess. The ratio of the energy density in the synchrotron photons. 
Us, to that in the magnetic field is. 



Us 



4 p-2 Ky] "a-j 



(45) 



B^'/Sk 3;r p - 3 0j//A//rj ' 

where AH is estimated as above, and the ratio is 4 for p = 4, 
7i = 10' and our adopted parameters. Thus, the self-Compton pro- 
cess is important, though Comptonization of blackbody radiation 
still dominates the electron losses, {Us + B-/Sn)/Uy^i - 0.2 [which 
reduces 7b, equation l ll9t . to =^ 10"] at rjB = 10"'. Note that since 
the electron distribution is determined by the observed 7-ray spec- 
trum (which is due to blackbody scattering), the above ratio is in- 
dependent of B. 

We assume that the synchrotron emission is isotropic in the jet 
frame. The Compton process is here mostly in the Thomson limit, 
and we treat it using a delta-function approximation, e = a'y^eg, 
where eg and e are the seed and scattered photon energy, respec- 
tively, and a' = I. This yields, 

.1/2 r™"<'/«)„s(eo), 



edhsc 
dedO. 



1/2 



— I deo s 



(46) 



(47) 



where iis is the density of the synchrotron photons, 

d«s d^is 4 

«s(fn) - 4:rr ^ cm \ 

deodn cV deadO. c&^HAHT^ 

and the integration limits in equation ( 146b account for the Thomson 
limit (assuring that 6 < 7) and 72 < 7 < 1. We have found that 
using a' = i reproduces better the exact Thomson-limit results for 
power law electrons (see Appendix iBt for 2 < p < 4 than a' = 4/3, 
corresponding to the average scattered energy. The resulting spec- 
tra, from both the jet and the counterjet, and taking into account the 
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relativistic transformation of equation M5\ . are shown in Figs. [3 a- 
b). We see tliat the self-Compton component may contribute to both 
the X-ray high-energy tail and the y-ray spectrum above 0. 1 GeV. 
The relative strength of this component is constrained by the orbital 
modulation. Intrinsically, the self-Compton component is not or- 
bitally modulated, which implies its contribution to the GeV range 
(with strong modulation) is weak. We find that the contribution of 
the self-Compton component to the range > 0. 1 TeV (not shown in 
Figs.[5^-b) is below the extrapolation of the FermifLAT power law. 
In the hard X-rays, we also see a relatively strong orbital modula- 
tion due to wind absorption jZdziarski et al.ll2012h . which appears 
to imply that the tail is not mainly due to this process. These con- 
straints on the relative amplitude of the self-Compton component 
also give an upper limit on B, as its increase would amplify both 
the synchrotron and self-Compton components by ~ 6'f +'"2 xfms^ 
we obtain i? < 10^ G within the y-ray emitting source. 

We also calculate the radial Thomson optical depth of the elec- 
trons, Tt = N^I(2n@jH^HTj), which is ^ 7 x 10"^ for our adopted 
parameters. 



6 DISCUSSION 

6.1 Electron acceleration 

It remains unclear what mechanism is responsible for acceleration 
of electrons producing high energy gamma-rays in jets. It is of- 
ten considered to be diffusive shock accel eration (PSA), which 
involves the first order Fermi pr ocess (e.g., iBlandford & Ostriked 
ll978l : lBlandford & EichleJigsV ). Initial studies of DSA scenarios 
were focused on explaining the origin of cosmic rays. The process 
was intensively explored by Monte Carlo simulations to take into 
acco unt different shock paramet ers and magnetic field structures 
(e.g. lNiemiec & Ostrowskj|2004l and references therein). Such sim- 
ulations fully confirmed the ability of DSA to produce ultrarela- 
tivistic cosmic rays, but acceleration of electrons (and positrons) 
was achieved only after assuming that they were preheated up to en- 
ergies corresponding with the momentum of thermal, shocked ions. 
Recent results obtained using particle in cell (PIC) simulations have 
shown that in collisionless relativistic electron-ion shocks a quasi- 
Maxwellian distribution of el ectrons is produced, with the aver- 
age energy of yim "Z ('Mi/'«e)7i dSironi &SDitkovskvll201ll, see also 
ISpitkovskvll200i) ." where % is the average Lorentz factor of ions. 
A quasi-Maxwellian distribution of preheated electrons with the 
temperature close to that of the io ns is also found in PIC simu - 
lations of non-relativistic shocks (.Riguelme & Spitkovskvll201 lb . 
This provides the required preheating, solving the above long- 
standing problem. 

There are still no available results of PIC simulations of mildly 
relativistic shocks. However, given the results mentioned above for 
relativistic and non-relativistic shocks, we can assume that a quasi- 
Maxwellian distribution of efficiently preheated electrons is pro- 
duced also in mildly relativistic shocks, which is likely to be the 
case in the Cyg X-3 jet. Noting that the relative contribution of 
the low-energy tail of a Maxwellian distribution to observed elec- 
tromagnetic spectra is small, we have approximated the electron 
injection spectrum to have a cut-off at an energy of yi ~ 7i„j. Then, 
our finding of yi ~ 300-10' is consistent with the above results, 
with this value being related to the m-Jm^ mass ratio. 

We note that low-energy breaks at y, » 1 are common in jets 
of AGNs, where the y are also often att r ibuted to the ion mass (e.g., 
IStawarz et all |2007|) . iGhisellini et"^ feOlOl. l201lh find that the 



electron distribution of blazars observed by Fermi commonly show 
a relatively steep injection above yi ~ 10^-10' (denoted in their pa- 
pers by yb), and a hard injection below it, for which the steady-state 
distribution is approximately compatible with that of our equation 
( I21l (. The low-energy cut-offs/breaks at yi ~ 10^-10' are observed 
not only in blazars but also in spectra of hot spots in radio -lobes 
telundell et alj2006l . [stawarz et alj2007l : lGodfrev et alj200 9) and 
used to argue for the presence of protons. Finally, we note that since 
the mass per ion in Cyg X-3 (a helium system) is ~3 times higher 
than for the cosmic abundances, the value of yi in it may be cor- 
respondingly higher than in comparable systems with abundances 
dominated by hydrogen. 

6.2 Caveats 

DCHIO obtained some ranges of the allowed parameters, but, given 
the complexity of the problem, we have just used the best-fit values 
(adjusted for the case of a steady jet) in this study. Furthermore, 
the best-fit parameters of DCHIO may be modified if the minimum 
of the X-ray folded light curve does not exactly correspond to the 
superior conjunction. This may happen if the wind is not symmetric 
with respect to the conjunctions, e.g., due to the wind lagging the 
binary rotation, and/or a formatio n of a Compton cloud ar ound the 
compact object, as in the model of lZdziarski et alj ( l2010al) . 

There is then a considerable uncertainty regarding the bi- 
nary parameters of Cyg X-3. We have adopted the binary pa- 
rameters for the black-hole case used by DCHIO, in particular 
the mass of the WR sta r of SO Mf?), which they assumed follow- 
ing Iszostek^^dzi^s^ ( 1200 Sl) . We note that such high mass ap- 
pears inconsistent with the mass-loss rate of ~ 10"''MQ/yr, es- 
timated for Cyg X-3 (e.g., Szostek & Zdziars ki 2008 and refer- 
ences therein). The mass vs. mass-loss rate relationship in WR 
stars is M ~ 10~^(M./Mo)"'M o/yr, with m 2-2.5 dLanged 
ll989l ; ISchaerer & Maedetlll992h . This implies M. ~ lOM p^, which 
also agrees with the results of Lomme n et al.l ( l2005h (though 
iHanson. Still & Fender 2000 favour a higher mass in the black-hole 
case). However, the value of M, affects only relatively slightly the 
orbital separation, oc (M, + M)'''. 

The stellar luminosity adopted by DCHIO, =^ 1.8 x lO^' erg 
s"', is close (within a factor of 2) to tha t predicted by the WR mass- 
luminosity relation tSchaerer & Maeder.. 1992 ). The chosen temper- 
ature, lO*" K, corresponds to the hydrostatic stellar surface rather 
than the photosphere of an isolated star, wit h the eff'ective tem- 
peratu re of the photosphere a few times lower JSchaerer & Maede^ 
1 19921) . However, the X-ray source in Cyg X-3 strongly ionizes the 
wind on the side of the jet, and thus the jet is likely to be exposed 
to radiation at the temperature close to the core one. On the other 
hand, Compton scattering of the stellar radiation by the wind will 
be substantial, which will increase the apparent size of the diluted 
blackbody source. 

In this study, we have used the Klein-Nishina cross section 
for calculating spectra, but still assumed the steady-state electron 
distribution is a power law. We note that such an approach is not 
fully self-consistent if the energy losses of the electrons by Comp- 
ton scattering dominate over the synchrotron an d adiabatic ones 
JZdziarski & Kroliklll993l: iModerski et al.ll2005[) . If the electrons 
are accelerated at a power-law rate, the electron energy losses re- 
duced in the Klein-Nishina regime cause the steady-state electron 
distribution to be no longer a power law. In fact, the reduced energy 
losses are largely compensated by the reduced Compton scattering 
emission, and the final spectrum is approximately a power law with 
the same index and normalization as that in the Thomson regime 
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jZdziarski & Kroliklll993l : lModerski et"ai]|2005h . Note that this ef- 
fect would need to be taken into account in calculating constraints 
on the maximum accelerated energy, 72. 

Finally, we have found that a model with a low value of the 
low-energy cut-off in the electron distribution, 71 < yb, and the 
break energy due to energy losses at yj, ^ 10' can also explain 
the y-ray spectrum of Cyg X-3. Although we cannot rule out this 
model, we consider it unlikely given our estimate of yb - 10^. 

6.3 The uniqueness of y-ray emission of Cyg X-3 

We briefly address the question why Cyg X-3 is, so far, the only 
accreting X-ray binary with confirmed high-energy y-ray emis- 
sion. We note that although a number of other X-ray binaries, 
e.g., LS I -l-61°303 or LS 5039 emit high-energy y-rays, that 
emission is, most likely, due to collision of their pulsar winds 
with stellar winds of their high-mass companions rather than due 
to accretion (e.g., Dubus '2006b'; 'Neron oy & Chernvakoval IiOOTI : 
IZdziarski. Neronov & Chernyakova 2010lJ! 

DCHIO noted that since Cyg X-3 has both a very high wind 
mass loss rate and a very small separation, it may be unique in 
forming a reconfinement shock in its inner jet. We note that, for 
the adopted Cyg X-3 parameters, this requires a relatively large ini- 
tial jet opening angle for a reconfinement shock to occur, 0j > 30°. 
This follows from equation (7) of DCHIO, which implies H > R 
for a smaller 0j, while H < R from the geometry of the system. 
Such a large initial opening angle may be for med in the jet forma- 
tion mechanism utilizing disc magnetic field (Blandford & Pavng 
Il982i) . We note that a si milar initial wide opening jet angle is seen 
in the radio galaxy M87 jSiretta. Junor & Livioll2002l) . 

A related unique feature of Cyg X-3 is the very large lumi- 
nosity of its companion (L, ~ 10^' erg s"') accompanying its 
very small separation. This results in a very high blackbody flux 
irradiating the y-ray emitting region, which then yields a strong 
Compton-scattering flux. If the orbital separation were > 10 times 
higher (as in Cyg X-1) or the stellar luminosity were much lower 
(as in low-mass X-ray binaries), synchrotron and self-Compton 
emission would dominate instead of blackbody up-scattering (see 
Section|53}. The presence of such synchrotron component may be 
searched for in those systems. The relative strength of the associ- 
ated self-Compton component in jets of those objects during dissi- 
pation events remains unknown (as it depends on the unknown jet 
parameters); we note it might produce observable y-ray emission. 



6.4 Relationstiip to radio emission 

After the relativistic electrons lose their energy at // ~ 10'^ cm, 
the jet continues to propagate for a large distance until another 
dissipation region forms. Since the stellar emission is very weak 
at that point, the main energy losses are synchrotron. The char- 
acteristic size of the sources of the resulting flaring radio emis- 
sion is ~ 10'^ c m, and the variability is on a day time scale (e.g., 
iMiller- Jones et al. 2004). In Fig. [T] we see that the average radio 
luminosity during the active periods is much lower than the y-ray 
luminosity. Thus, the observed radio emission is energetically al- 
lowed to be emitted by the jet downstream the y-ray dissipation re- 
gion, with its power reduced by about an order of magnitude (Sec- 
tion |5.1t with respect to the upstream jet. 

If the jet experiences a dissipation episode of the kind stud- 
ied here, its velocity in the radio-emitting region should be sim- 
ilar to that in the dissipation region (estimated from the orbital 



modulation), jSj ~ 0.5. Inde ed, such a velocity has be en estimated 
from the proper motion by iMiller- Jones et aL I ( I2OO6I) , and a sim- 
ilar ySj ^ 0.6 was estimated from fitting a pr ecession model to 
a radio image obtained during a flaring state jMiller -Jones et al.l 
,2004.) . Similar estimates have been obtained from a number of 
other radio observations of Cyg X-3, e.g., ( 3. ^ 0.5 inferred by 
iMartf. Paredes & Peracaul j j200l l). except for Mioduszewski et all 
1 200 ih , who estimated j8j > 0.8. We note that our estimate of the jet 
velocity before the dissipation region is compatible w ith that, jSi p ~ 
0.9. It is then possible that the measurement of Mioduszewski et al.l 
( l200lh was done during a radio flaring episode during which the 
dissipation region producing y-rays was not formed, and the jet 
propagated to large distances with /3j.o. Occurrence of radio flaring 
episodes without formation of a prior dissipation region at scales 
<K 1 0'^ cm can also expla in a radio flare occurring before a y-ray 
flare jWiUiams etai]|201ll) . 

Recently, mm radio flares lasting a fraction of a day and oc- 
currin g on intermediate size sca les, ~ 10'^ cm, have been discov- 
ered i TTsuboi et~al.ll20ld I2OI2I) . This size scale is derived from 
the flare rise time scales of several minutes. The second flare of 
iTsuboi et ai] ( l2012h . observed at 43 GHz and 86 GHz, took place 
on MJD 54972.9 during a period quiescent in both radio emission 
at 15 GHz and in y-rays (FLC09). Thus, this is an example of a 
radio flare without a prior jet energy dissipation on the orbital size 
scale. The electron power required was found > 3 X (Z)/7 kpc)- 10" 
erg s"' (neglecting relativistic corrections), similar to our estimate 
of the electron power. The magnetic field was estimated at ~ 10 G 
(assuming equipartition), an order of magnitude below our upper 
limit for the y-ray emitting region. 



7 CONCLUSIONS 

We have obtained the average X-ray spectrum and the radio flux 
emitted during the y-ray active epochs (Section O. We have then 
calculated spectra from Compton scattering of a photon beam into 
a given direction by isotropic relativistic electrons with a power- 
law distribution with a low-energy cut-off. Simple analytical for- 
mulae have been obtained both using the Klein-Nishina cross sec- 
tion (Section[3) and in the Thomson limit (AppendixlAt. 

We have applied our results to scattering of stellar blackbody 
radiation by relativistic electrons in the jet of Cyg X-3 (Section|4]l, 
using the model of DCHIO, fitted to the observed modulation of 
y-rays (FLC09). We have found a low-energy break at yi ~ 300- 
10' in the distribution of the accelerated electrons is required by the 
observational data in order not to overproduce the observed X-ray 
emission. We find Compton cooling to be eflicient, which implies 
the power-law index of the acceleration process of ^ 2.5-3, rather 
typical to astrophysical acceleration sites. The low-energy electron 
break found by us is in agreement with recent shock acceleration 
models, in which it is related to the ion/electron mass ratio. Also, 
the obtained value of the break Lorentz factor is similar to those 
typically found in AGN jets (see Section [6T| l. 

We have calculated the jet kinetic power to be ~ 10'^ erg s"' 
assuming equipartition between the energy supplied to the electrons 
and ions, which represents a firm lower limit. The estimated power 
is comparable to the radiative power of this source. Assuming this 
equipartition, the bulk Lorentz factor of the jet before the dissipa- 
tion region is ~ 2.5. Most of the power supplied to the electrons is 
radiated (Sections [STTlfS^ . 

We have found that the magnetic field strength is constrained 
to be below the equipartition with the electron energy density by 
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a factor of < a few times 10"^. At the upper limit (5 ~ 10^ G), 
the synchrotron emission from the y-ray emitting region still gives 
rise to a relatively strong IR flux, which measurement simultaneous 
with that of y-rays would provide an estimate of the magnetic field 
in this part of the jet (Section [5.3t . The predicted synchrotron flux 
is at the level at the IR flux measured during past radio flares of 
Cyg X-3. 
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APPENDIX A: THE THOMSON LIMIT 

In the Thomson hmit, the flux per electron and per solid angle be- 
comes, 
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where, in order to approximately constrain the resulting spectrum 
to the energy range satisfying the Thomson limit, we have imposed 
a sharp cut-off at the range boundary. This, in particular, assures 
e < 7. Then, the condition of r < 1 yield the electron-integrated 
rate given by equation l|5j but with the lower and upper limits of 
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respectively. These limits, the constraint of 6 < 260x72, and a 
Thomson limit condition of e < l/(2x6o) yield, 
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When 6 < 7J/72, there is no range for /i - /4, and when e> y\ly'i, 
there is no range for fi - fi- When 72 — > 00, only the f\ - and 
,/2 - .A ranges exist. When e > 72, the spectrum is null. The above 
formulae should be applied for e » 60 only. The emitted spectrum 
in the part dominated by fi is a power law with an energy index of 
a = {p- l)/2. 

Figs.[2]and[3]shows compares the spectra obtained using equa- 
tions ( IA4b and Q for blackbody irradiation, respectively, with the 
corresponding ones obtained using the Klein-Nishina cross section. 
When 6/(27^x) « <6o> <c min[e/(27jx), l/(2xe)], where (6o> is 
the characteristic energy of the seed photons, the integral over /i 
in equation (O dominates, and we can set its limits from zero to 
infinity. 

In the case of diluted blackbody seed photons, equation l ll2b . 
and for e/(27^x) <c kT/{m^c^D,) « min[6/(27^x), l/(2xe)], the 
scattered flux in a given direction per dimensional energy in the jet 
frame becomes. 



edh 
dedn 



2^3o-TC~h-'K(em^c')^(kT)- 



ll+4p + p- /l+p\ /5+p 
5+p \ 2 



(A9) 



where ^ is the Riemann function. This flux was obtained by 
DCHIO, see their equations (1), (3). [Their formula is in the ob- 
server's frame assuming the 2)| transformation, see Section|4l it is 
47t times larger than ours, which results from their definition of K 
as corresponding to dA'/d7dn rather than dA'/d7, and the power of 
2 in their equation (3) is misprinted as p + 5/2, while is should be 
(p + 5)/2.] 



APPENDIX B: ISOTROPIC SCATTERING 

We note that integrating spectra over the scattering angle gives ei- 
ther the spectrum for a photon beam integrated over all directions of 
the scattered photon or the spectrum from scattering on an isotropic 
seed photon distribution. Also, the case of isotropic seed photons 
and isotropic electrons is equivalent to an electron beam when scat- 
tered photons are integrated over all directions. 

Therefore, integrating equation l lAlt over cos 1? in the range 
corresponding to r < 1 and multiplying by 2n (corresponding to in- 
tegration over the azimuth), we o btain the isotropic rate of equation 
(2.42) of iBlumenthal & Gouldl ( Il970l) . Integrating of equation 
l lA4t over x and multiplying by 2;r gives power-law electron rates 
integrated over all directions, which also correspond to emission 
at any direction in the case of isotropic seed photons. In particu- 
lar, assuming that we are far below both the maximum emitted e 
by electrons with 72 and the boundary of the Thomson limit, we 
obtain. 



de 87T " \/r(e,eo,p,7iX e < 4eo7f, 



fr = 2n f Mx = 2^*"n 
Jo 



e > 4eo7?; 
? < 4eo7f, 

\l+4p + p- 



(Bl) 



(B2) 



Jo 



-<fF/(2eoyj) /^2 

: 2;r I fidx + 2n 

Jc/acojh 



{I + p){3 + p)H5 + p) 

/2dx = -^x (B3) 



867t 2e^y] 



fo(l+p) 4(3 + p) 



1 +21n 



4eo7? 



4(3 + p)^ ehS+p) 



These rates can also be obtained by integrating th e isotropic rate 
of equation (2.42) of IBlumenthal & Gouldl imd) over the elec- 
tron distr ibution Q. The rate o f equat ion iB2i is given by equation 
(2.64) o f IBlumenthal & Gould JT970l) . We can then integrate over 
a photon distribution. 
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Integrating the rate of equation ( IA9I ) over x from to 2 and 
multiplyin g by 2;r gives the iso tropic rate of equations (2.65-2.66) 
of Blumen thal & Gouldl ( Il970h ). except that their rate is given for 
the seed photon density inside a blackbody field, with no(eo) in 
units of erg"' cm"', rather for the photon flux, as in our equa- 
tion I ll2t . That rate can also be obtained from equation dB4t for 
6/(47,) ^ kT ID, (i.e., neglecting the integral over ff°). Integrat- 
ing the Klein-Nishina rate of equation l|6ll over the scattering angle 
le ads to the isotropic Klein-Nis hina rate of equations (25) and (A9) 
of lAharonian & Atovarj ( ll98lh . 



© 2012 RAS, MNRAS OOO.fTirB] 



